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ABSTRACT 

We report the discovery and characterization of four transiting exoplanets by the HATNet survey. 

The planet HAT-P-50b has a mass of 1.35 Mj and radius of 1.29 i?j, and orbits a bright {V = II.8 mag) 

M = 1.27 Mq, R = I.7Oi?0 star every P = 3.1220days. The planet HAT-P-51b has a mass of 0.31 Mj 
and radius of 1.29Pj, and orbits a U = 13.4mag, M = 0.98 Mq, R = 1.04Pq star with a period 
of P = 4.2180days. The planet HAT-P-52b has a mass of 0.82 Mj and radius of l.OlPj, and orbits 
a U = 14.1 mag, M = 0.89M©, R = 0.89P© star with a period of P = 2.7536days. The planet 
HAT-P-53b has a mass of 1.48Mj and radius of 1.32Pj, and orbits a U = 13.7mag, M = 1.09M©, 

P = 1.21 P© star with a period of P = 1.9616 days. All four planets are consistent with having 
circular orbits and have masses and radii measured to better than 10% precision. The low stellar 
jitter and favorable Pp/P* ratio for HAT-P-51 make it a promising target for measuring the Rossiter- 
McLaughlin effect for a Saturn-mass planet. 

Subject headings: planetary systems — stars: individual ( HAT-P-50, GSC 0787-00340, HAT-P- 
51, GSC 2296-00637, HAT-P-52, GSC 1793-01136, HAT-P-53, GSC 2813-01266 ) 
techniques: spectroscopic, photometric 
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1. INTRODUCTION 


Transiting exoplanets (TEPs) are important objects 
for studying the physical properties of planets outside 
the solar system. By combining time-series photome¬ 
try of a transit with time-series radial velocity (RV) ob¬ 
servations of the star spanning the planetary orbit, it 
is possible to accurately measure the mass and radius 
of a transiting planet relative to those of the host star. 
Leveraging stellar evolution models to estimate the stel¬ 
lar mass and radius given observable parameters such as 
the effective temperature, metallicity and bulk density 
of the star, then allows the physical mass and radius of 
the planet, as well as its orbital separation, to be deter¬ 
mined. Other properties of the sys tem such as the or - 
bital eccentricity and obliquity le.g. lQueloz et al.ll200fll l. 
and properties of the planetary atmosphere (e.g. emis¬ 
sion or transmission s pectra) may also be access ible for 
transiting planets fe.g. iCharbonneau et al.l[2002D . Moti¬ 
vated by the wealth of physical information that may be 
measured for these objects, there has been a significant 
effort over the past 15 years to discover and character¬ 
ize many TEPs. The aim of this effort is to explore the 
diversity of exoplanets, and to identify statistically ro¬ 
bust relations between their physical parameters, which 
in tu rn inform theories of planet formation and evolution 
('e.g.lGu illot et al']l2006HBnrrows et al.ll200' 7l: iBekv et al.l 

I2OI1I: iLaughlin et al.ll2011l: lEnoch et al.ll201^ 


Largely thanks to the ultra-high-precision photometric 
time-series observations from the NASA Kepler mission, 
we now know of over 4000 high-quality ca ndidate tran¬ 
siting exoplanets fe.g. iMullallv et al]l2015D . Some 51 of 
the Kepler candidates have been confirmed through mea¬ 
suring the RV orbital wobble of their host stars, while a 
further 845 have masses estimated through transit time 
variations, or have been statistically validated as being 
ven^ unlikely to be anything other than transiting plan- 


iry u 

etqij. The majority of the candidates from Kepler are, 
however, too small and/or orbiting stars that are too 
faint to allow their masses and orbital eccentricities to 
be determined using existing spectroscopic facilities. Eor 
most of these planets, all we can determine at present are 
their radii, orbital periods, and a constraint on their ec- 
centricities using the so-c alled photo-eccentric effect (e.g. 
IDawson fc JohnsoDll2012T) . 

Most of the TEPs with spectroscopically determined 
masses have been discovered by wid e-held ground-base d 
transit surv eys such as HAT Net (iBakos e t al.l I20 04D. 
HAT South (iBakos et al.l l201p). W ASP (iPollacco et al 


2006 ), XO (iMcCullough et al.ll2005l). Tr ES ( Alonso et al. 


200 I . and KETT (jPenner et al.m^ . among others. 


These surveys cover a greater area of the sky than has 
been surveyed so far by Kepler or its successor mission 
K2, and have thereby monitored more bright stars which 
may host TEPs amenable to conhrmation spectroscopy. 
In this paper we present the discovery and characteriza¬ 
tion of four new transiting short-period gas-giant planets 
by the HATNet survey. 

The HATNet survey, which began operations in 2004, 
has to date searched 17% of the 47r steradian celestial 
sphere for planets. A total of 5.5 million stars have 
been observed. The stars have from 2400 to 21000 high- 
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cadence photometric observations (5th and 95th per¬ 
centiles; the median is 7200) spanning a few months 
to several years. The point-to-point RMS precision of 
the observations ranges from ~ 3 mmag for stars with 
r ^ 9 to ~ 2% for stars with r ^ 13.3 (depending on 
sky conditions and the density of stars in the field being 
observed). Based on these observations we have iden¬ 
tified ^ 2000 candidate TEPs, the majority of which 
are false positives. The stars are generally bright (the 
median magnitude of the candidates is V = 12.7mag) 
so that it has been possible to carry out spectroscopic 
and/or photometric follow-up observations for the ma¬ 
jority of these objects. Based on this follow-up, 1468 
candidates have been rejected as false positives (the tran¬ 
sit signal is probably real, but not due to a planet), 189 
have been rejected as false alarms (the identified tran¬ 
sit signal was not real), while more than 50 confirmed 
and well-characterized planets (including those presented 
here) have been announced. Some ^ 350 candidates are 
currently active. 

The four new planets announced in this paper have 
properties that are typical of short-period gas-giant plan¬ 
ets. While they do not, in themselves, reveal new proper¬ 
ties of exoplanets, they will contribute to our statistical 
understanding of planetary systems in the Galaxy. 

In the next section we describe the observations used to 
confirm the new TEPs. In Section [3] we discuss the anal¬ 
ysis carried out to rule out false positive blend scenarios 
and determine the physical parameters of the planetary 
systems. We place these planets into context with the 
other known transiting planets in Section IT] 

2. OBSERVATIONS 

The discovery of all four transiting planet systems fol- 
lowed the gen e ral ob serva tional procedures de scribed by 
iLatham et al.l (|2009D and IBakos et al.l (I2010D . Here we 
summarize the observations of each system, and our 
methods for reducing the raw data to scientifically in¬ 
teresting measurements. 

2.1. Photometric detection 

The four TEPs presented here were initially identi¬ 
fied as candidate TEPs based on observations made 
with the HATNet wide-held photometric instruments 
(IBakos et al.l 1200^ . This network consists of six iden¬ 
tical fully-automated instruments, with four at Pred 
Lawrence Whipple Observatory (PLWO) in AZ, and two 
on the roof of the Submillimeter Array Hangar Build¬ 
ing at Mauna Kea Observatory (MKO) in HI. The light¬ 
gathering elements of each instrument include an 11 cm 
diameter telephoto lens, a Sloan r hlter, and a 4Kx4K 
front-side-illuminated CCD camera. Observations made 
in 2007 and early 2008 were done using a Cousins R hl¬ 
ter. The instruments have a held of view of 10? 6 x 10? 6 
and a pixel scale of 9"pixel“^ at the center of an image. 
Observations are fully automated with the typical pro¬ 
cedure being to continuously monitor a given held while 
it is above 30° elevation taking exposures of 180 s (prior 
to 2010 December an exposure time of 300s was used). 
The helds have been dehned by tiling the sky into 838 
7° X 7° pointings. Because each tile is smaller than the 
held of view, there is overlap between neighboring helds 
and a given source may be observed in multiple (up to 
four) helds. 
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Table [T] lists the HATNet observations which con¬ 
tributed to the discovery of each system. All four ob¬ 
jects were observed using multiple HAT instruments, and 
three of the four objects are in overlapping fields. In some 
cases the observations date back to 2007, and may span 
as many as 3.5 years. HAT-P-51, in particular, has been 
observed extensively with HATNet, having more than 
27,000 individual photometric measurements. 

The raw HATNet images were re duced to syst e matic- 
noise-filtered light curves following iBakos et al.l (|2010[1 
and making use of a perture an d image subtraction pho¬ 
tometry tools from iPall (j2009f l. The filtering includes 
decorrelating the individual light curves against various 
instrumental parameters (we refer to this procedure as 
External Parameter Decorrelation, or EPD) including 
the image position of the source, the sub-pixel position, 
the background flux, the local scatter in the background 
flux, and the shape of the PSF. Following EPD we make 
use o f the Trend Filtering Algorithm (TEA: iKovacs et aP 
in non-reconstructive mode. The data for each 
HATNet field were reduced independently, with EPD ap¬ 
plied separately to each instrument, and TEA applied 
globally to all observations from a given field (with an 
option to perform a complete TEA filtering, using data 
from all telescopes and all fields). 

Light curves were searched for periodic box-shaped 
transi ts using the Box-fi tting Least Squares algorithm 
fBLS: li^vacs et ;iI1l2nn2[ l. Candidates were selected us¬ 
ing a variety of automated cuts (e.g. on the S/N, dif¬ 
ferences in depth between even and odd transits, among 
others) and a final by-eye inspection. Figured] shows the 
phase-folded, trend-filtered light curves from HATNet for 
the four newly discovered planetary systems. 

We used BLS to search the residual light curves for 
additional planetary transits, but did not detect any ad¬ 
ditional signals. We als o calculated the Discrete Fourier 
Transform (PET, see iDee min a EM and using the 
method of iKurtd 119851 for a fast recursive evaluation of 
the trigonometric functions) for each of the light curves, 
after subtracting the best-fit transit models, to search 
for any continuous periodic variations. Such variations 
may be due to the rotation of spotted stars, for exam¬ 
ple. For HAT-P-50, -52 and -53 we can rule out signals 
in the frequency range 0 to 50 d“^ with an amplitude 
above O.Gmmag, l.Smmag and 1.3mmag, respectively. 
For HAT-P-51 we also do not find a significant Fourier 
component. Curiously, the highest peak in the frequency 
spectrum is within 1.3% of the first harmonic of the or¬ 
bital frequency. We do not have a physical explanation 
of this near coincidence, if it is a real signal, but we can 
exclude the possibility of tidal distortion due to the well- 
demonstrated sub-stellar nature of the companion (see 
Section [32]). It may perhaps be a signature of stellar ac¬ 
tivity. After subtracting this low amplitude (l.Smmag) 
component, the next highest peak in the frequency spec¬ 
trum has an amplitude of l.Ommag. 

2 . 2 . Spectroscopic Observations 

Follow-up spectroscopic observations were carried out 
using six different facilities. The aim of these observa¬ 
tions was to aid in ruling out false positives, determine 
the atmospheric parameters of the host stars, and to 
confirm the planets by measuring the RV orbital vari¬ 
ations induced by the transiting planets. The facili¬ 


ties used for each system are summarized in Table [5] 
and include t he Tillinghas t Reflector Echelle Spectro¬ 
graph ITRES: lFuresd[200§l on the 1.5 m Tillinghast Re¬ 
flector at FLWO; the Astrophysical Research Consor - 
tium Echelle Spectrometer lAR.CES: iWang et al.ll2003h 
on the ARC 3.5 m telescope at Apache Point Observa¬ 
tory (APO) in New Mexico; the Fibre-fed Echelle Spec¬ 
trograph (FIES) at the 2.5 m Nordic Optical Telescop e 
(NOT) at La Palma, Spain (jDiuDvik fc Andersenll201(lll : 
the S OPHIE Spectrograp h on the 1.93 m t elescope at 
OHP (|Bouchv et al.ll200^ in France; HIRES (|Vogt et al.l 
1199411 on the Keck-I telescope in Hawaii together with the 
I 2 abso rption cell; and the H igh-Dispersion Spectrograph 
HDS: iNo giic li e t al.l 1200^ with the I 2 absorption cell 
Kambe et al. I2002fl on the Subaru telescope in Hawaii. 

The TRES observations were used for reconnaissance 
(i.e. ruling out false positives with lower S/N spectra) 
for HAT-P-51, HAT-P-52 and HAT-P-53. For HAT-P-50 
they were used both for reconnaissance and for measuring 
the orbital variation due to the planet. The raw echelle 
images were reduced to extracted spectra and analyzed 
to me asure RVs and stellar a tmospheric parameters fol¬ 
lowing [Buchh^e^E^ (|lQi3). Observations of standard 
stars were made during each observing run and are used 
to correct the velocities from each run to the lAU system. 
Because these corrections are known for TRES, we adopt 
the TRES measurements for the systemic 7 velocity of 
each object listed in Table [4] The uncertainty on the 
absolute calibration is ^ 0 . 1 kms“^ and is dominated by 
the uncertainty in the absolute velocities of the standard 
stars. 

The ARCES observations of HAT-P-51 and HAT-P- 
53 were used exclusively for reconnaissance (based on 
observations of standard stars the RV precision of this 
instrument is limited to ^ 500ms“^). Observations 
were reduced to wavelength calibrated spectra using the 
ECHELLE package in IRAll^. For the wavelength calibra¬ 
tion we made use of ThAr lamp spectra obtained before 
or after each science exposure, and with the same point¬ 
ing as the science exposure. Each spectrum was ana¬ 
lyzed to measure the RV of the star, its surface gravity, 
effective temperature, projected equatorial rotation ve¬ 
locity, and meta llicity using the Stella r Parameter Clas¬ 
sification (SPC: iBiichhave et al.ll20lll procedure, which 
cross-correlates the observed spectrum against a set of 
synthetic template spectra. 

The single FIES spectrum obtained for HAT-P-51 was 
used for reconnaissance, and wa s reduced and analyzed 
following iBuchhave et al.l (l201Cl[l . 

SOPHIE observations of HAT-P-51 were collected in 
high-efficiency mode with the aim of confirming the 
planet by measuring the RV orbital wobble of its host 
star. The SOP HIE observat i ons w ere reduced and ana¬ 
lyzed following iBoisse et al.l (1201311 . Based on these ob¬ 
servations we determined that HAT-P-51b is a Saturn- 
mass planet, and that the 40ms“^ precision of the 
SOPHIE observations for this object was insufficient to 
accurately determine the planetary mass. The precision 
in this case was limited due to significant contamination 

IRAF is distributed by the National Optical Astronomy Ob¬ 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 1. Phase-folded unbinned HATNet light curves for HAT-P-50 (upper left), HAT-P-51 (upper right), HAT-P-52 (lower left), and 
HAT-P-53 (lower right). In each case we show two panels. The top panel shows the full light curve, while the bottom panel shows the light 
curve zoomed-in on the transit. The solid lines show the model fits to the light curves. The dark filled circles in the bottom panels show 
the light curves binned in phase with a bin size of 0.002. 


from scattered moon light, for uncontaminated spectra 
significantly higher precision may be obtained from the 
same S/N. We do not include these data in the analysis 
of HAT-P-51. 

HDS observations were collected for HAT-P-50 and 
HAT-P-51 in order to confirm these TEP systems and 
characterize the planetary orbits. The observations were 
extracted and reduced to relati ve RVs in the solar sys - 
tem barycentric frame following iSato et al.l (I2002L f2(Tl^ . 
while spe ctral line bis e ctor s pans (BSs) were computed 
following [Torres et al.l (1200711 . 

HIRES observations were collected for HAT-P-51, 
HAT-P-52 and HAT-P-53. The observations have an RV 
precision of 5-10 ms“^ and are used here to character¬ 
ize the orbital variations and to determine the stellar 
atmospheric parameters. The data were re duced to rela¬ 
tive R Vs in the barycentric frame following iButler et al.l 
(|1996ll . Spectra l-line bisecto r spans (BSs) were com¬ 
puted following [Torres et al.l (12007), and S activity 
indice s were calculated following [Isaacson fc Fischeil 
(|2010D These were transf ormed to log^Q values fol¬ 
lowing |NoyeiIeOn (|1984D . 

Based on the reconnaissance TRES, FIES and ARCES 
observations we find that none of the four targets shows 
evidence of being a composite system. All have radial ve¬ 
locity variations below lkms“^, and all are dwarf stars. 
The effective temperatures, projected rotation velocities, 
and surface gravities estimated from these spectra are 


consistent with the higher precision values presented in 
Table 01 

The high-precision RV measurements for all objects are 
seen to vary in phase with the transit ephemerides. These 
are shown in Figure 0J In this same figure we also show 
the phased BS measurements, which in all cases are con¬ 
sistent with no variation in phase with the ephemerides. 
The data are listed in Table [Hat the end of the paper. 

2.3. Photometric follow-up observations 

Additional time-series photometric measurements were 
obtained for all four of the systems using Keplercam on 
the FLWO 1.2 m telescope. These observations were car¬ 
ried out during the planetary transits to aid in ruling out 
blended eclipsing binary false positive scenarios, and to 
refine the light curve parameters (i.e. the orbital period, 
the planet to star radius ratio, the impact parameter and 
the transit duration). For HAT-P-50 we also obtained 
follow-up photometry with the CCD imager on the Byrne 
Observatory at Sedgwick (BOS) 0.8 m telescope, located 
at Sedgwick Reserve in Santa Ynez Valley, CA, and oper¬ 
ated by the Las C umbres Observatory Global Telescope 
institute (LCOGT: [Brown et al.[[^13l) . The events mon¬ 
itored with each instrument, together with the number 
of images obtained, the cadence, filter used and photo¬ 
metric precision are listed in Table [T] 

We applied standard CCD calibration procedures to 
the Keplercam and BOS images and then reduced these 
to light curves using the aperture photometry methods 
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Table 1 

Summary of photometric observations 


Instrument/Field^ 

Date(s) 

# Images 

Cadence'^ 

(sec) 

Filter 

Precision*^ 

(mmag) 

HAT-P-50 






HAT-10/G316 

2008 Nov-2009 May 

3214 

352 

Sloan r 

7.5 

HAT-5/G364 

2009 May 

21 

411 

Sloan r 

10.6 

HAT-9/G364 

2008 Dec-2009 May 

3159 

352 

Sloan r 

7.5 

BOS 

2012 Feb 15 

105 

149 

Sloan i 

2.1 

Keplercam 

2012 Feb 18 

443 

54 

Sloan i 

1.2 

BOS 

2012 Feb 21 

81 

140 

Sloan i 

2.5 

BOS 

2012 Apr 08 

61 

143 

Sloan i 

1.6 

Keplercam 

2012 Nov 28 

462 

44 

Sloan i 

1.8 

Keplercam 

2012 Dec 23 

277 

45 

Sloan i 

2.3 

Keplercam 

2013 Jan 14 

427 

45 

Sloan i 

1.4 

Keplercam 

2013 Jan 17 

380 

45 

Sloan i 

1.6 

HAT-P-51 






HAT-6/G164 

2007 Sep-2008 Feb 

3652 

349 

Cousins R 

30.3 

HAT-9/G164 

2007 Sep-2008 Feb 

2767 

349 

Cousins R 

25.9 

HAT-10/G165 

2010 Sep-2011 Jan 

4215 

230 

Sloan r 

24.3 

HAT-5/G165 

2010 Nov-2011 Feb 

4142 

354 

Sloan r 

24.1 

HAT-8/G165 

2010 Nov-2011 Feb 

2240 

238 

Sloan r 

23.6 

HAT-6/G209 

2010 Nov-2011 Feb 

3794 

351 

Sloan r 

18.4 

HAT-9/G209 

2010 Nov-2011 Feb 

2151 

352 

Sloan r 

18.0 

HAT-7/G210 

2010 Nov-2011 Jan 

4047 

229 

Sloan r 

19.1 

Keplercam 

2011 Oct 21 

88 

134 

Sloan i 

1.9 

Keplercam 

2012 Jan 05 

92 

133 

Sloan i 

2.7 

Keplercam 

2012 Oct 05 

171 

134 

Sloan i 

2.2 

Keplercam 

2012 Oct 26 

137 

134 

Sloan i 

2.6 

Keplercam 

2012 Nov 12 

111 

134 

Sloan i 

3.2 

HAT-P-52 






HAT-5/G212 

2010 Sep-2010 Nov 

2270 

347 

Sloan r 

19.5 

HAT-8/G212 

2010 Aug-2010 Nov 

5999 

232 

Sloan r 

22.4 

Keplercam 

2010 Dec 23 

101 

134 

Sloan i 

2.0 

Keplercam 

2011 Sep 05 

90 

133 

Sloan i 

2.7 

Keplercam 

2011 Sep 27 

188 

134 

Sloan i 

2.3 

Keplercam 

2011 Nov 21 

82 

133 

Sloan i 

2.5 

Keplercam 

2012 Jan 07 

64 

194 

Sloan i 

3.0 

HAT-P-53 






HAT-6/G164 

2007 Sep-2008 Feb 

3653 

349 

Cousins R 

26.4 

HAT-9/G164 

2007 Sep-2008 Feb 

2764 

349 

Cousins R 

24.5 

HAT-10/G165 

2010 Sep-2011 Jan 

4234 

230 

Sloan r 

19.3 

HAT-5/G165 

2010 Nov-2011 Feb 

4134 

354 

Sloan r 

19.4 

HAT-8/G165 

2010 Nov-2011 Feb 

2240 

238 

Sloan r 

20.4 

Keplercam 

2011 Oct 19 

158 

134 

Sloan i 

1.9 

Keplercam 

2011 Oct 27 

381 

73 

Sloan i 

2.5 


^ For HATNet data we list the HAT station and field name from which the observations are taken. 

HAT-5, -6, -7, and -10 are located at FLWO in Arizona, while HAT-8 and -9 are located at MKO 
in Hawaii. Each field corresponds to one of 838 fixed pointings used to cover the full 47r celestial 
sphere. All data from a given HATNet field are reduced together, while detrending through External 
Parameter Decorrelation (EPD) is done independently for each unique field-|-station combination. 

^ The mode time between consecutive images rounded to the nearest second. Due to weather, the 
day—night cycle, guiding and focus corrections, and other factors, the cadence is only approximately 
uniform over short timescales. 

^ The RMS of the residuals from the best-fit model. 

described bv IBakos et al.l (I2010D . In doing this we made the four objects are available in machine-readable form 
use of the stellar centroid positions measured directly in Table [3l 
from a set of registered and stacked frames rather than 
relying on catalog positions for astrometry as done in 
iBakos et alJ (1201011 . All sources in the images, save the 
target TEP system, were used in performing the ensem¬ 
ble magnitude calibration. We corrected for additional 
systematic trends in the data by including the EPD and 
TEA noise filtering models in the fitting mentioned in 
Sectioning The resulting trend filtered light curves for 
HAT-P-50 through HAT-P-53 are shown in Figures [3l- 
1 respectively. All photometric measurements made for 
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Table 2 

Summary of spectroscopy observations 


Instrument 


Date(s) 

# Spec. 

Res. 

AA/A/1000 

S/N Range 

RV Precision^ 
(ms~^) 

HAT-P-50 

TRES 

2010 

Dec-2012 Feb 

5 

44 

24.8-36.1 

116.0 

FIES 

2012 

Mar 13-17 

5 

67 

31.0-69.9 

25.0 

HDS 

2012 

Feb 7 

3 

60 

271-283 


HDS+I2 

2012 

Feb-2012 Sep 

20 

60 

84-166 

23.0 

HAT-P-51 

FIES 

2011 

Aug 4 

1 

46 

27.4 


ARCES 

2011 

Sep 19 

1 

31.5 

20.6 


TRES 

2011 

Sep 21 

1 

44 

20.9 


SOPHIE 

2011 

Dec 4-12 

4 

39 

23-28 

37 

HIRES 

2011 

Oct-Nov 

2 

55 

83-94 


HIRES+I2 

2011 

Oct-2012 Feb 

6 

55 

59-80 

5.4 

HDS 

2012 

Feb 9 

4 

60 

52-56 


HDS+I2 

2012 

Feb 7-10 

20 

60 

26-53 

9.2 

HAT-P-52 

TRES 

2010 

Dec-2011 Jan 

2 

44 

19.1-20.4 

300 

HIRES 

2011 

Oct 19 

1 

55 

66 


HIRES+I2 

2011 

Feb-2012 Jul 

7 

55 

26-59 

7.5 

HAT-P-53 

TRES 

2011 

Sep 18-19 

2 

44 

30.6-30.7 

80 

ARCES 

2011 

Sep 19-20 

2 

31.5 

19.8-20.1 

380 

HIRES 

2011 

Nov 14 

1 

55 

90 


HIRES+I2 

2011 

Nov-2012 Feb 

6 

55 

62-79 

11.0 


^ The signal-to-noise ratio per resolution element near 5180 A. 

^ The RMS of the RV residuals from the best-fit orbit, or the RMS of the RVs for reconnaissance 
observations. We do not give an estimate for template spectra (listed as HIRES or HDS without 12 
included), or for cases where only a single spectrum was obtained with a given instrument. 
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HAT-P-50 


HAT-P-51 




0.4 0.6 
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Figure 2. Phased high-precision RV measurements for HAT-P-50 (upper left), HAT-P-51 (upper right), HAT-P-52 (lower left), and 
HAT-P-53 (lower right). In each case we show three panels. The top panel shows the phased measurements together with our best-fit 
model (see Table Ell for each system. Zero-phase corresponds to the time of mid-transit. The center-of-mass velocity has been subtracted. 
The second panel shows the velocity O — C residuals from the best fit. The error bars include the jitter terms listed in Table [5] added in 
quadrature to the formal errors for each instrument. The third panel shows the bisector spans (BS), with the mean value subtracted. The 
symbols used for each instrument are indicated in the top panel for each planet. Note the different vertical scales of the panels. 
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HAT-P-50 
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Figure 3. Left: Unbinned transit light curves for HAT-P-50, acquired with the CCD imager on the B OS 0 .8 m telescope, and Keplercam 
on the FLWO 1.2 m telescope. The light curves have been EPD and TFA processed, as described in fj 13.31 The dates of the events and 
instrumen ts us ed are indicated. Curves after the first are displaced vertically for clarity. Our best fit from the global modeling described 
in Section 13.31 is shown by the solid lines. Right: Residuals from the fits in the same order as the left panel. The error bars represent the 
photon and background shot noise, plus the readout noise. 
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HAT-P-51 



Figure 4. Similar to Figure [3] here we show the follow-up light curves for HAT-P-51. All light curves were obtained with Keplercam on 
the FLWO 1.2 m telescope. 
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HAT-P-52 



Figure 5. Similar to Figure |3] here we show the follow-up light curves for HAT-P-52. All light curves were obtained with Keplercam on 
the FLWO 1.2 m telescope. 
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Figure 6. Similar to Figure |3] here we show the follow-up light curves for HAT-P-53. All light curves were obtained with Keplercam on 
the FLWO 1.2 m telescope. 
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3. ANALYSIS 

3.1. Properties of the parent star 

The stellar atmospheric parameters that we adopted 
for the analysis, including the effective temperature Teff*, 
the surface gravity logg*, the metallicity [Fe/H] and the 
projected equatorial rotation velocity usini, were deter¬ 
mined for each system using SPC. For HAT-P-50 we ap¬ 
plied this to the TRES and FIES spectra (applying to 
the individual spectra and adopting the average parame¬ 
ter values) while for the other three systems we used the 
Keck/HIRES l2-free template sp ectra. _ 

We used the Yonsei-Yale (Y2; lYi et al.l 1200111 theoret¬ 
ical stellar models to determine physical parameters of 
the stars, such as their masses, radii, luminosities and 
ages, based on the measured atmospheric parameters to¬ 
gether with the bulk stellar densities p* determined from 
our modelling of the light curves and RV measurements 
(Section [231) • We generated a chain of Teff*, [Fe/H] and 
p* values for each object, where the p* values are taken 
from the output of the MCMC procedure used to fit 
the light curves and RVs, while we assume uncorrelated 
Gaussian distributions for Tgff* and [Fe/H]. For each 
value in the chain we interpolate the Y2 models to find 
a combination of M*, age and [Fe/H] which matches the 
three input parameters (we assume solar-scaled abun¬ 
dances without a-element enhancement). Combinations 
of Teff*, [Fe/H] and p* that do not match to a stellar 
model are rejected. In doing this we also reject the cor¬ 
responding link in the LC-I-RV MCMC chain so that the 
final planetary parameters are restricted to regions of 
parameter space allowed by the stellar evolution models. 
The stellar models also provide other parameters such 
as i?* and L* for a given M*, age and [Fe/H] combina¬ 
tion. The result is a posterior chain of stellar parameters 
for each star. We use the chains to calculate the me¬ 
dian and 68.3% confidence interval for each of the stellar 
parameters. These are listed in Table [4| We compare 
the measured Tgff* and p* values for each system to the 
model isochrones in Figure [7| 

For HAT-P-50 and HAT-P-52 we found that the me¬ 
dian logp* values determined from this procedure dif¬ 
fered significantly from the values estimated from the 
spectra. For these stars we carried out a second itera¬ 
tion of SPC fixing logp* to the values determined from 
the stellar evolution modelling. We then performed a 
second iteration of the LC-I-RV modelling, with revised 
limb darkening parameters, followed by a second itera¬ 
tion of the stellar evolution modelling. The log p* values 
had converged after this iteration. For HAT-P-51 and 
HAT-P-53 a second iteration of SPC was not needed. 

Distances are determined for each system by compar¬ 
ing the measured broad-band photometry listed in Ta¬ 
ble [4] to the magnitudes predicted in each filter by the 
models. We allow for extinction assumi ng a Ry = 3.1 
extinction law from iCardelli et ahi (|1989t) . 

Based on this modelling we find that HAT-P-50 has a 
mass of 1.273^0.115^10, a radius of 1.698 ± 0.071 i?©, an 
age of 3.37lo 27 Gyr, and is at a distance of 497 ± 21 pc. 
HAT-P-51 has a mass of 0.976 ± 0.028M©, a radius 
of 1.041^0029-^0) of 8-2 i 1-7 Gyr, and is at 

a distance of 470 ± 16 pc. HAT-P-52 has a mass of 
O.887±O.O27M0, a radius of 0.893±0.047i?©, an age of 


9.4 ±4.1 Gyr, and is at a distance of 385 ±21 pc. Finally, 
HAT-P-53 has a mass of 1.093 ± 0.043M0, a radius of 
1.209lo;og2 i?0, an age of 4.67^083 and is at a dis¬ 
tance of 719 ± 43 pc. 

For HAT-P-51, -52 and -53 we used the Keck/HIRES 
spectra to determine median logig Rhk activity indices. 
We find that all three stars are inactive in the Ca H HK 
region, consistent with their slow rotation and lack of 
photometric variability. 

3.2. Excluding blend scenarios 

In order to excl ude blend scenarios w e carried out an 
analysis following ([Hartman et al.l[2012fl . We attempt to 
fit the available photometry (light curves, and catalog 
broad-band magnitudes calibrated to standard systems) 
for each object using a combination of three stars (two 
eclipsing, with a third diluting the eclipse signal) with 
properties taken from stellar evolution models. 

For HAT-P-50, HAT-P-51 and HAT-P-53 we find that 
a model consisting of a planet transiting an isolated star 
provides a better (lower y^) fit to the data than any of 
the blend models tested. For HAT-P-50 the best-fit blend 
model is excluded with 1.5 ct confidence, while for both 
HAT-P-51 and HAT-P-53 it is excluded with 2a con¬ 
fidence. We also simulated cross-correlation functions, 
RVs and BS measurements for the blend models tested, 
and found that any model that could plausibly fit the 
photometry for these systems (i.e. provides a fit that is no 
more than 5cr worse than the single star±planet model) 
would be easily identified as a composite stellar system 
based on the spectroscopy. We therefore conclude that 
all three of these objects are transiting planet systems. 

For HAT-P-52 we similarly find that the planet±star 
model provides a better fit to the data than any blend 
model tested, however the best-fit blend model differs 
by only 0.5cr from the planet±star model. We also find 
that there is a range of models consisting of a blend be¬ 
tween a bright foreground star, and a background stel¬ 
lar eclipsing binary that is between 0.35 mag and 4 mag 
further in distance modulus than the foreground star, 
which cannot be ruled out based on the photometry or 
BS spans. For these models the simulated BS variations 
have a scatter that is below the 43ms“^ scatter in the 
Keck/HIRES data, if we allow for a difference in the 7 
velocities of the foreground star and background binary. 
We find, however, that similar to the case of HAT-P- 
49 ([Biervla et al.l[20T^ . the expected form of the RV 
variations in these blends is significantly different from 
the observed sinusoidal variation, even though the over¬ 
all amplitude of the variations is comparable (Fig. [8]). 
We conclude that although the photometry and BS mea¬ 
surements for HAT-P-52 can be fit by a blended stellar 
eclipsing binary model, the RV observations cannot be. 

While we can rule out the possibility that any of these 
objects is a blended stellar eclipsing binary system, we 
cannot rule out the possibility that one or more of these 
transiting planet systems also has a stellar companion. 
For HAT-P-50, we find that models including a faint 
companion with M > 0.6 Mq provide a slightly worse 
fit to the data than models without a companion. The 
difference in is small, however, and we can only rule 
out a stellar companion with M > 1.2 Mq at greater than 
3(7 confidence. For HAT-P-51, models with a companion 
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Table 3 

Light curve data for HAT-P-50-HAT-P-53. 


Object 

BJO*^ 

(2,400,000-1-) 

Mag'^ 

o-Mag Mag(orig)‘^ 

Filter 

Instrument 

HAT-P-50 

54863.83699 

0.00928 

0.00393 

r 

HATNet 

HAT-P-50 

54935.64419 

-0.00034 

0.00517 

r 

HATNet 

HAT-P-50 

54838.86248 

0.00069 

0.00443 

r 

HATNet 

HAT-P-50 

54888.81494 

-0.00915 

0.00442 

r 

HATNet 

HAT-P-50 

54792.03291 

0.00668 

0.00450 

r 

HATNet 

HAT-P-50 

54910.67022 

0.00400 

0.00418 

r 

HATNet 

HAT-P-50 

54863.84110 

0.00001 

0.00388 

r 

HATNet 

HAT-P-50 

54935.64826 

0.00737 

0.00652 

r 

HATNet 

HAT-P-50 

54838.86690 

0.00070 

0.00433 

r 

HATNet 

HAT-P-50 

54888.81903 

0.00426 

0.00441 

r 

HATNet 


Note. — This table is available in a machine-readable form in the online journal. A portion 
is shown here for guidance regarding its form and content. 

^ Either HAT-P-50, HAT-P-51, HAT-P-52, or HAT-P-53. 

^ Barycentric Julian Date is computed directly from the UTC time without correction for leap 
seconds. 

^ The out-of-transit level has been subtracted. These magnitudes have been subjected to the 
EPD and TEA procedures, carried out simultaneously with the transit fit. 

Raw magnitude values without application of the EPD and TEA procedures. These are provided 
only for the follow-up observations. For HATNet, the transits are only detectable after applying 
the noise filtering methods. 
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HAT-P-52 



Effective temperature [K] 


HAT-P-53 



Effective temperature [K] 


Figure 7. Model isochrones from lYi et al.l Il2nnii^ for the measured metallicities of HAT-P-50 (upper left), HAT-P-51 (upper right), 
HAT-P-52 (lower left) and HAT-P-53 (lower right). In each case we show models for several different ages, with younger models being 
on the left and older ones on the right. For HAT-P-50, HAT-P-51 and HAT-P-52 we show ages of 0.2 Gyr and 1.0 to 14.0 Gyr in 1.0 Gyr 
increments. For HAT-P-53 we show ages of 0.2 Gyr and 1.0 to 9.0 Gyr in 1.0 Gyr increments. The adopted values of T^ff* and p* are shown 
together with their Itr and 2a confidence ellipsoids. The initial values of Tefj* and p* from the first SPG and light curve analyses are 
represented with a triangle for HAT-P-50 and HAT-P-52. Note the logarithmic vertical axes, and the different scales used in each panel. 


having M > 0.5 M0have a slightly worse fit to the data, 
but we can only rule out companions with M > 0.95 Mq 
at greater than 3 ct confidence. For HAT-P-52, compan¬ 
ions with M > 0.35 Mq provide a slightly worse fit to the 
data, but all companions up to the mass of HAT-P-52 are 


permitted to within 3 ct. For HAT-P-53 companions with 
M > 0.7 Mq provide a slightly worse fit, but all compan¬ 
ions up to the mass of HAT-P-52 are permitted to within 
3cr. 
























14 


Hartman et al. 


Table 4 

Stellar parameters for HAT-P-50-HAT-P-53 



HAT-P-50 

HAT-P-51 

HAT-P-52 

HAT-P-53 


Parameter 

Value 

Value 

Value 

Value 

Source 

Astrometric properties and cross-identifications 




2MASS-ID. 

07521521+1208218 

01241564+3248387 

02505320+2901206 

01272906+3858053 


GSC-ID. 

GSC 0787-00340 

GSC 2296-00637 

GSC 1793-01136 

GSC 2813-01266 


R.A. (J2000). 

07*'52“15.20s 

0l''24™ 15.66s 

02‘'50™53.20s 

0l‘'27™29.05s 

2MASS 

Dec. (J2000). 

+12°08'21.9" 

+32°48'38.8" 

+29°01'20.6" 

+38°58'05.3" 

2MASS 

Mr.A. (masyr-l) 

10.20 + 0.80 

-9.8 + 1.4 

12.5 + 2.1 

-1.1 ± 1.8 

UCAC4 

MDec. (masyr-l) 

-4.6 ±1.6 

-16.8 + 2.1 

-24.7 + 2.6 

3.0+ 2.2 

UCAC4 

Spectroscopic properties 






TefiMK). 

6280 ± 49 

5449 ± 50 

5131 + 50 

5956 ± 50 

SPC'^ 

[Fe/H]. 

-0.180 ±0.080 

0.270 + 0.080 

0.280 ± 0.080 

0.000 ± 0.080 

SPC 

V sin i (km s ^) ... 

8.90 + 0.50 

1.70 + 0.50 

0.60 + 0.50 

4.10 + 0.50 

SPC 

Vmac (kms“l)- 

1.0 

1.0 

1.0 

1.0 

Assumed 

Vmic (kms-1) .... 

2.0 

2.0 

2.0 

2.0 

Assumed 

7rv (kms-1). 

6.25 + 0.10 

-27.56 + 0.10 

61.50 + 0.10 

-16.99 + 0.10 

TRES 

logic ^HK ^ 


-5.057 + 0.050 

-5.154 + 0.089 

-4.919 ± 0.042 

HIRES 

Photometric properties 






B (mag). 

12.282 + 0.050 

14.261 ± 0.070 

15.183 + 0.050 


APASS,TASS'^ 

V (mag). 

11.762 + 0.030 

13.440 ± 0.040 

14.068 + 0.020 

13.73 + 0.18 

APASS,TASS'= 

Ic (mag). 

11.194 + 0.052 

12.67 + 0.12 

13.02 + 0.17 

13.13 ± 0.12 

TASS 

9 (mag). 

11.973 + 0.050 

13.839 ± 0.050 

14.631 + 0.060 


APASS 

r (mag). 

11.650 + 0.030 

13.194 + 0.030 

13.677 + 0.080 


APASS 

* (mag). 

11.550 + 0.020 

12.998 ± 0.040 

13.441 + 0.090 


APASS 

J (mag). 

10.816 + 0.021 

12.039 ± 0.022 

12.195 + 0.022 

12.468 ± 0.023 

2MASS 

H (mag). 

10.545 + 0.020 

11.645 ± 0.023 

11.745 + 0.022 

12.202 ± 0.026 

2MASS 

Ks (mag). 

10.500 + 0.018 

11.614 ± 0.020 

11.621 + 0.021 

12.100 + 0.019 

2MASS 

Derived properties 






(Mq) . 

1 Q'7Q”i"0-049 

0.976 ± 0.028 

0.887 + 0.027 

1.093 + 0.043 

YY+a/R*+SPC 

R* (Re) . 

1.698 + 0.071 

029 

0.893 ± 0.047 

Qg2 

YY+a/R*+SPC 

log 9* (cgs). 

4.072 + 0.029 

4.392 ± 0.027 

4.483 + 0.051 

4.310 + 0.043 

YY+a/R*+SPC 

P* (gcm-3). 

0.357 + 0.037 

1 OQo+U.lUU 
i.zzo_o 135 

1.75 + 0.29 

0.87 + 0.13 

YY+a/R*+SPC 

M (Lq) . 

4.01 + 0.38 

0.859 ± 0.070 

0.496 + 0.060 

1 65+0-24 
J--do_q 

YY+a/R*+SPC 

Mv (mag). 

3.27 + 0.11 

5.055 ± 0.095 

5.73 + 0.13 

4.27 + 0.13 

YY+a/R*+SPC 

Mk (magjESO) . . 

2.052 + 0.093 

3.268 ± 0.071 

3.69 + 0.12 

2.83 + 0.12 

YY+a/R*+SPC 

Age (Gyr). 

o 07+1-44 
-0.27 

8.2 ± 1.7 

9.4+ 4.1 

4 67+^+^ 
^•O'_o.83 

YY+a/R*+SPC 

Ay (mag). 

0 011+0-056 

U.Uil_o_oii 

u.uiz_o 012 

0.412 + 0.052 

0.21 + 0.14 

YY+a/R*+SPC 

Distance (pc). 

497 + 21 

470 ± 16 

385 ± 21 

719 ± 43 

YY+a/R*+SPC 


^ SPC — “Stellar Parameter Classification” routine for the analysis of high-resolution spectra IIBuchhave et ani2012ll . applied to the TRES 
and FIES spectra of HAT-P-50, and to the Keck/HIRES l 2 -free template spectra of HAT-P-51, HAT-P-52 and HAT-P-53. These parameters 
rely primarily on SPC, but have a small dependence also on the iterative analysis incorporating the isochrone search and global modeling of 
the data, as described in the text. 

^ The median of the logj^o ^hk values measured from the individual Keck/HIRES spectra for each target. The uncertainty is the standard 
error on the median. 

From APASS DR6 for HAT-P-50, HAT-P-51 and HAT-P-52 as listed in the UCAC 4 catalog iZacharias et aI1l2012fl . From TASS Mark IV 
jDroege et al.|[2006ll for HAT-P-53. 

^ YY-ho/R*-hSPC = Based on the YY isochrones er et alJl2n0ll . a/Ri, as a luminosity indicator, and the SPC results. 


3.3. Global modeling of the data 

We modeled the HATNet photometry, the follow-up 
photometry, and the high-precision RV measurements us - 
ing the procedur e descr ibed in detail bv iPal et al.l (l2008f l 
andIBak os et al. (1201011 with modifications described by 
iHartman et al.l ( 201211 . This procedure makes use of the 
differential evol ution Markov Chain Monte Carlo (DE M- 
CMC) method (|ter Braakll2006l : lEastman et al.l[2013ll to 
explore the fitness landscape and produce posterior pa¬ 
rameter distributions. We allowed for RV jitter which 
we varied as a free parameter in the fit for each planet. 
We adopted independent jitters for each instrument as 
the methods for estimating the “formal” errors differ be¬ 
tween reduction methods and instruments. For HIRES 
we made use of an empirical prior on the jitter as dis¬ 


cussed in IHartman et all (|2014f) . while for the other in¬ 
struments we used a Jeffreys prior (i.e. the prior proba¬ 
bility for parameter cr is oc 1/cr). We fixed the limb dark - 
ening coefficients using the tabulation by iClare 3 (l200l 
and the stellar atmospheric parameters given in Tabled 
The resulting parameters for each system are listed 
in Table [5l We find that HAT-P-50b is a hot Jupiter 
with a mass of 1.350 ± 0.073 Mj and radius of 1.288 ± 
0.064 i?j, HAT-P-51b is a hot Saturn with a mass of 
0.309 ± 0.018Mj and radius of 1.293 ± 0.054i?j, while 
HAT-P-52b and HAT-P-53b are hot Jupiters with masses 
of 0.818 ± 0.029 Mj and 1.484 ± 0.056 Mj, and radii of 
1.009 ± 0.072 i?j and 1.318 ± 0.091 i?j, respectively. We 
fit all systems both allowing the eccentricity to vary and 
fixing it to zero. We find that all four systems are con¬ 
sistent with no eccentricity (the 95% confidence upper 
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Phase 

Figure 8 . Phase-folded Keck/HIRES RVs for HAT-P-52 com¬ 
pared to the best-fit planetary orbit (dashed curve) and the best- 
fit RVs from a blended stellar eclipsing binary model (solid curve). 
This same blend model has BS variations that are below the ob¬ 
served scatter, and fits the photometric data, however it does not 
reproduce the form of the RV variation, which is well described by 
a transiting planet. The blend-model fit shown here consists of a 
0.88 Mq foreground star blended with a 0.84 Af0-^O.13 Mq eclips¬ 
ing binary that has a distance modulus 0.35 mag greater than the 
foreground star, and which has a 7 velocity that is 70kms“^ differ¬ 
ent from that of the foreground star. Reducing the difference in 7 
velocities creates a more symmetric RV variation, but also results 
in the RV variation going to 0ms“^ at both phases 0.25 and 0.75. 
All other blend models simulated have RV variations that provide 
even poorer fits to the observations. 

limits on the eccentricity when it is allowed to vary are 
e < 0.115, < 0.123, < 0.047, and < 0.134 for HAT- 
P-50b through HAT-P-53b, respectively). We therefore 
adopted the parameters for a fixed circular orbit in all 
cases. 

4. DISCUSSION 

In this paper we presented the discovery and charac¬ 
terization of four transiting exoplanets from the HAT- 
Net survey, including three hot Jupiters (HAT-P-50b, 
HAT-P-52b and HAT-P-53b) and a hot Saturn (HAT- 
P-51b). All four planets have masses and radii deter¬ 
mined to better than 10% precision. The mass uncer¬ 
tainties are 5.4%, 5.8%, 3.5%, and 3.8% for HAT-P-50b 
through HAT-P-53b, respectively, while the respective 
radius uncertainties are 5.0%, 4.2%, 7.1%, and 6.9%. 
The stars HAT-P-50, -51, and -53 also have fairly pre¬ 
cise isochrone-based age determinations (uncertainty less 
than 2 Gyr) thanks to their favorable position within the 
Tes-p* plane (Fig.[71). 

In Figure m we show the location of these planets on a 
mass-radius diagram, comparing them to the full sam¬ 
ple of confirmed TEPs with 0.1 < Mp < lOMj. The 
new planets all fall within the range of values already 
seen by other planets, with HAT-P-51b falling near the 
upper envelope of the distribution of points in the mass- 
radius diagram, and HAT-P-52b falling near the lower 
envelope. We also show the location of each planet on 
a Teq-radius diagram. Again we find that the planets 
all follow the well-established trends. While not atypi¬ 
cal compared to other known exoplanets, these objects 
contribute to the growing sample of well-characterized 
planets which may be used to explore the population of 
planets in the Galaxy through statistical methods. 

In terms of potential for additional follow-up observa¬ 
tions, we conclude that it should be feasible to measure 
the Rossiter-McLaughlin effect for HAT-P-50b, HAT-P- 


51b and HAT-P-53b using Subaru/HDS or Keck/HIRES. 
For HAT-P-50b the expected amplitude of the R-M ef- 
fect is 42m s~^ for an aligned orbit (using eq. 40 in 
IwiTiTllarinh . For HAT-P-51b the expected amplitude 
is 27ms“^, and for HAT-P-53b it is 48ms“^. The Sub¬ 
aru/HDS velocity residuals for HAT-P-50 have an RMS 
of 23ms“^, with a median exposure time of 10 minutes. 
Assuming 20 such exposures are obtained over the course 
of a single transit, it should be possible to measure the 
R-M amplitude to a precision of 5a (based on fitting 
models to simulated observations). For HAT-P-51, the 
Keck/HIRES RVs have a residual RMS of 5.4 ms“^, and 
a median exposure time of 25 minutes. Seven of these ex¬ 
posures could be collected over a single transit, allowing 
a 9a detection of the R-M amplitude. For HAT-P-53, 
the Keck/HIRES RVs have a residual of 11 ms“^ and an 
exposure time of 25 minutes. For this system it should 
be possible to collect 6 similar exposures during a tran¬ 
sit, and measure the R-M amplitude with Scr confidence. 
For HAT-P-52b the R-M amplitude is only < 7ms“^ 
(limited by the very slow rotation), and we would not 
expect to detect it in a single transit with better than 2a 
confidence. 

The conclusion that the R-M effect should be easier to 
detect for both HAT-P-51 and HAT-P-53 than for HAT- 
P-50, despite both stars being significantly fainter than 
HAT-P-50, and despite both stars having a lower nsini, 
may be counter intuitive. The RV observations for HAT- 
P-51 and HAT-P-53 are both significantly higher preci¬ 
sion than those for HAT-P-50, even though the HAT-P- 
50 observations have higher S/N. Some of the difference 
may be due to the different instruments (Subaru/HDS 
for HAT-P-50 vs. Keck/HIRES for HAT-P-51 and HAT- 
P-53). However, slower rotation and cooler surface tem¬ 
peratures are also factors which tend to improve the RV 
precision. In this respect we expect HAT-P-51 to have 
higher precision than HAT-P-53 at fixed S/N, and HAT- 
P-53 to have higher precision than HAT-P-50 at fixed 
S/N, which is what we see. 

Measuring the R-M effect for HAT-P-51b may be 
of particular interest due to its small mass. HAT-P- 
11b and Kepler-63b are the only planets smaller than 
HAT- P-51b for which t his ef fect has been measured to 
date (jWinn et al.l 120101 and iSanchis-Oieda et Ml 120131 
the obliquity has also been measur ed for the Kepler-30 
syste m by star-spot crossings, see iSanchis-Oieda et ^ 

[20T11 . 

While the R-M effect should be detectable for HAT-P- 
50b, HAT-P-51b and HAT-P-53b, due to the relatively 
small value of Rp/ i?* for HAT-P-50b, and the faintness of 
the other targets, none of the new planets are particularly 
well-suited for atmospheric characterization. 

HATNet operations have been funded by NASA grants 
NNG04GN74G, NNX08AF23G, and NNX13AJ15G. 
Follow-up of HATNet targets has been partially sup¬ 
ported through NSF grant AST-1108686. G.A.B, Z.C. 
and K.P. acknowledge partial support from NASA grant 
NNX09AB29G. J.H. acknowledges support from NASA 
grant NNX14AF87G. K.P. acknowledges support from 
NASA grant NNX13AQ62G. G.T. acknowledges partial 
support from NASA grant NNX14AB83G. We acknowl¬ 
edge partial support also from the Kepler Mission un- 
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Table 5 



Orbital and planetary parameters for HAT-P-50b 

-HAT-P-53b 


Parameter 

HAT-P-50b 

Value 

HAT-P-51b 

Value 

HAT-P-52b 

Value 

HAT-P-53b 

Value 

Light curve parameters 

P (days) . 

Ta (BJD) . 

3.1220109 ±0.0000065 
... 2456285.90993 ± 0.00036 

4.2180278 ± 0.0000059 
2456194.12204 ± 0.00040 

2.7535953 ± 0.0000094 
2455852.10326 ± 0.00041 

1.9616241 ± 0.0000039 
2455829.44781 ± 0.00044 


Ti 4 (days) 

Ti 2 = T 34 (days) . 

a/Ri, . 

C/R* . 

Rp/R* . 

{,2 . 

b = a cos i / it* . 

i (deg) . 

Limb-darkening coefficients 

Cl, i (linear term) . 

C 2 ,i (quadratic term) . 


0.1531 ±0.0011 
0.0176 ±0.0013 
5.68 ±0.19 
14.710 ±0.077 
0.0782 ± 0.0012 

U.OaO_Q ggg 
U.DZa_g g 4 j^ 

83.65 ±0.57 


0.1965 

0.3570 


0.1403 ±0.0016 
0.0170 ±0.0012 
10.48i°j« 
16.21 ±0.12 
0.1278 ±0.0020 

U.U/ I _0,g52 
n oyy-t-0.085 

88.48 ± 0.57 


0.3348 

0.2989 


0.1003 ± 0.0017 
0.0130 ± 0.0018 
8.89 ± 0.49 
22.88 ± 0.27 
0.1161 ± 0.0027 
n 91 Q-t"0.096 

U.410_g ggg 

0 461 +0 094 

U.4Di_g 422 
87.02 ± 0.86 


0.3908 

0.2628 


0.1164 ±0.0017 
0.0135 ±0.0016 
5.61 ±0.28 
19.42 ±0.17 
0.1120 ±0.0019 
0.142t°;°®9 
0.38l°;“ 
86.2 ± 1.5 


0.2387 

0.3447 


RV parameters 

K (ms-l) . 

e-l . 

RV jitter HIRES (ms"!) .. 

RV jitter HDS (ms“^) . 

RV jitter TRES (ms“^) . 

RV jitter FIES (ms"!) . 


161.3 ±5.6 
< 0.115 

0.7 ±7.1 
68 ±38 
0.00 ±0.91 


39.5 ±2.2 
< 0.123 
4.3 ±1.2 
0.0 ±1.4 


128.4 ±3.8 
< 0.047 
5.2 ± 1.9 


226.8 ±6.0 
< 0.134 
9.3 ±3.2 


Secondary eclipse parameters 

Ts (BJD) . 

Ts,i4 (days) . 

Ta,i 2 (days) . 

Planetary parameters 

Mp{Mj) . 

Rp (Rj) . 

C{Mp,Rp) f . 

Pp (gcm-3) . 

log 9p (cgs) . 

a (AU) . 

Teq (K) . 

08 . 

logio('l^) (cgs) ** . 


2456287.47093 ± 0.00036 
0.1531 ±0.0011 
0.0176 ±0.0013 


1.350 ±0.073 
1.288 ±0.064 
0.48 

0.78 ±0.11 
3.302 ±0.038 

U.Uq09U_g gg44g 

1862 ± 34 
0.0751 ± 0.0044 
9.433 ± 0.032 


2456196.23106 ± 0.00040 
0.1403 ±0.0016 
0.0170 ±0.0012 


0.309 ±0.018 
1.293 ±0.054 
-0.03 

0.178 ±0.024 

Z.DDl_g gg 4 

0.05069 ± 0.00049 
1192 ± 21 
0.0247 ±0.0018 
8.659 ±0.031 


2455853.48006 ± 0.00041 
0.1003 ± 0.0017 
0.0130 ± 0.0018 


0.818 ±0.029 
1.009 ± 0.072 
-0.15 
0.98 ± 0.21 
3.296 ± 0.065 
0.03694 ± 0.00038 
1218 ±37 
0.0673 ± 0.0050 
8.696 ± 0.052 


2455830.42863 ± 0.00044 
0.1164 ±0.0017 
0.0135 ±0.0016 


1.484 ±0.056 
1.318 ±0.091 
0.34 

0.80 ±0.15 
3.325 ± 0.055 
0.03159 ± 0.00042 
1778 ± 48 
0.0649 ± 0.0046 
9.354 ± 0.046 


^ Times are in Barycentric Julian Date calculated directly from UTC without correction for leap seconds. Tc". Reference epoch of mid transit that 
minimizes the correlation with the orbital period. T 14 : total transit duration, time between first to last contact; T 12 — T 34 : ingress/egress time, time 
between first and second, or third and fourth contact. 

^ Reciprocal of the half duration of the transit used as a jump parameter in our MCMC analysis in place of a/R*. It is related to a/R* by the expression 
C/R* — u/R*(27r(l -t- e sin a;))/(R-\/l — h^Vl — e^) dBakos et al.ll2010fl . 

Values for a quadratic law, adopted from the tabulations bv ICIar'^ II2004I1 according to the spectroscopic (SPC) parameters listed in Table HI 
As discussed in Sect ion 13.31 the adopted parameters for all four systems are determined assuming circular orbits. We also list the 95% confidence upper 
limit on the eccentricity determined when ^/ecosuj and y/esinuj are allowed to vary in the fit. 

® Term added in quadrature to the formal RV uncertainties for each instrument. This is treated as a free parameter in the fitting routine. For HIRES we 
include an empirical prior constraint following IHartman et al.l II2014I1 . 

^ Correlation coefficient between the planetary mass Mp and radius Rp. 

8 The Safronov number is given by 0 — ^ (Vesc/Vorb)^ — (n/Rp)(Afp/Af*) fsee lHansen &€ Barmanll2007[l . 

® Incoming flux per unit surface area, averaged over the orbit. 


der NASA Cooperative Agreement NCC2-1390 (D.W.L., 
PI). Data presented in this paper are based on observa¬ 
tions obtained at the HAT station at the Submillimeter 
Array of SAG, and the HAT station at the Fred Lawrence 
Whipple Observatory of SAG. The authors wish to ac¬ 
knowledge the very significant cultural role and reverence 
that the summit of Manna Kea has always had within the 
indigenous Hawaiian community. We are most fortunate 
to have the opportunity to conduct observations from 
this mountain. This research has made use of Keck tele¬ 
scope time granted through NO AO (program A245Hr) 
and NASA (N154Hr, NlSOHr). This research was made 


possible through the use of the AAVSO Photometric All- 
Sky Survey (APASS), funded by the Robert Martin Ay¬ 
ers Sciences Fund. 
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Figure 9. Left: Mass-radius diagram for transiting planets with 0.1 Mj < M < 10Mj. The four planets discovered here are indicated. 
The two filled gray squares show Saturn and Jupiter. The parameters for other transiting planets are compiled from the literature 
(c.f. http://www.exoplanets.org). Right: planet radius vs. estimated equilibrium temperature (assuming zero albedo and complete 
redistribution of heat) for the same sample of planets as shown at left. Symbols are assigned colors based on the planetary masses. 
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Table 6 

Relative radial velocities, bisector spans, and activity index measurements for 
HAT-P-50-HAT-P-53. 


BJD 


O'RV*’ 

BS 

<^BS 


Phase 

Instrument 

(2,454,000+) 

(ms“^) 

(ms“^) 

(ms (ms 




HAT-P-50 

1903.95158 

169.24 

31.82 




0.656 

TRES 

1911.98432 

-149.26 

31.82 




0.229 

TRES 

1941.89144 

-5.12 

65.07 




0.808 

TRES 

1957.72714 

81.84 

63.24 




0.881 

TRES 

1958.68285 

-19.27 

59.81 




0.187 

TRES 

1960.70621 

-26.06 

73.39 




0.835 

TRES 

1964.90560 

-161.79 

16.35 

-44.94 

22.34 


0.180 

HDS 

1964.91336 

-158.48 

15.07 

-40.90 

21.84 


0.182 

HDS 

1964.92112 

-164.86 

14.23 

-54.18 

20.02 


0.185 

HDS 

1964.93269 



-62.87 

21.17 


0.189 

HDS 

1964.94740 



-62.65 

21.82 


0.193 

HDS 

1964.96211 



-61.14 

22.32 


0.198 

HDS 

1965.06347 

-184.08 

15.46 

-110.01 

24.62 


0.231 

HDS 

1965.07123 

-196.36 

14.93 

-153.60 

27.34 


0.233 

HDS 

1966.86139 

143.00 

14.72 

-32.02 

20.00 


0.807 

HDS 

1966.87609 

137.51 

13.33 

-19.54 

19.70 


0.811 

HDS 

1967.86542 

-96.04 

21.28 

-7.35 

21.69 


0.128 

HDS 

1967.87318 

-93.64 

18.64 

-15.27 

20.80 


0.131 

HDS 

1967.88093 

-87.90 

16.20 

0.26 

18.99 


0.133 

HDS 

1967.88867 

-108.44 

17.37 

-9.93 

20.59 


0.135 

HDS 

2000.39479 

78.60 

23.25 




0.547 

FIES 

2001.40519 

121.33 

23.25 




0.871 

FIES 

2002.47018 

-191.77 

30.26 




0.212 

FIES 

2003.40044 

4.13 

54.05 




0.510 

FIES 

2004.38851 

120.16 

31.15 




0.827 

FIES 

2191.13792 

121.23 

22.53 

45.30 

19.61 


0.644 

HDS 

2191.14223 

129.48 

26.61 

129.85 

21.78 


0.645 

HDS 

2191.14652 

184.86 

29.91 

496.30 

30.08 


0.647 

HDS 

2192.12442 

48.71 

17.29 

7.33 

16.97 


0.960 

HDS 

2192.12871 

45.28 

23.29 

-3.33 

20.31 


0.961 

HDS 

2192.13300 

28.57 

21.63 

9.00 

17.56 


0.963 

HDS 

2193.11572 

-151.84 

16.13 

-7.12 

19.16 


0.277 

HDS 

2193.12002 

-121.63 

23.24 

-5.25 

17.21 


0.279 

HDS 

2193.12433 

-127.80 

20.99 

2.08 

16.32 


0.280 

HDS 

HAT-P-51 

1853.80027 

-42.75 

2.53 

-16.75 

5.54 

0.1550 

0.317 

HIRES 

1853.81798 



-2.36 

6.06 

0.1300 

0.322 

HIRES 

1879.85898 



10.56 

6.81 

0.1630 

0.495 

HIRES 

1879.87608 

7.25 

2.68 

8.32 

6.08 

0.1570 

0.499 

HIRES 

1880.97599 

42.01 

2.48 

9.54 

5.32 

0.0769 

0.760 

HIRES 

1903.95906 

-35.87 

2.63 

-42.11 

21.39 

0.1550 

0.209 

HIRES 

1944.73569 

24.93 

2.55 

-9.88 

7.15 

0.1170 

0.876 

HIRES 

1960.81411 

31.81 

3.23 

5.14 

7.02 

0.1600 

0.688 

HIRES 

1964.71582 

44.69 

13.17 

-0.33 

26.24 


0.613 

HDS 

1964.72359 

27.51 

11.37 

-23.77 

23.80 


0.615 

HDS 

1964.73475 

21.63 

11.55 

7.03 

29.96 


0.617 

HDS 

1964.74946 

25.14 

9.83 

3.26 

34.15 


0.621 

HDS 

1964.76417 

33.95 

9.59 

0.65 

35.27 


0.624 

HDS 

1965.71533 

32.14 

15.07 

-2.89 

24.12 


0.850 

HDS 

1965.73005 

29.21 

12.47 

9.32 

36.71 


0.853 

HDS 

1965.74477 

33.24 

11.95 

-14.30 

27.64 


0.857 

HDS 

1965.75948 

34.67 

14.11 

-16.82 

27.55 


0.860 

HDS 

1965.77418 

13.59 

12.75 

21.67 

22.10 


0.864 

HDS 

1965.78890 

15.87 

21.41 

23.88 

25.83 


0.866 

HDS 

1966.72018 



-4.21 

31.52 


0.088 

HDS 

1966.73489 



9.59 

32.96 


0.092 

HDS 

1966.74959 



4.25 

32.50 


0.095 

HDS 

1966.76429 



5.05 

34.50 


0.099 

HDS 

1966.77907 

-12.89 

12.16 

-0.33 

34.37 


0.102 

HDS 

1966.79377 

-33.08 

11.65 

27.17 

29.63 


0.106 

HDS 

1966.80848 

-13.58 

15.24 

5.11 

13.87 


0.109 

HDS 

1967.72003 

-24.46 

14.33 

-2.79 

22.98 


0.325 

HDS 

1967.73475 

-48.37 

13.90 

-12.59 

23.40 


0.329 

HDS 

1967.74946 

-42.32 

13.39 

-11.16 

21.89 


0.332 

HDS 

1967.76417 

-34.03 

12.75 

0.39 

20.57 


0.336 

HDS 

1967.77888 

-28.16 

13.62 

3.97 

17.43 


0.339 

HDS 

1967.79359 

-28.23 

14.33 

-32.17 

31.51 


0.343 

HDS 

HAT-P-52 
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Table 6 — Continued 


BJD 

RV==' 

o'RV*’ 

BS 

•^BS 

S‘= 

Phase 

Instrument 

(2,454,000+) 

(ms-i) 

(ms-l) 

(ms (ms 





1611.81789 

134.02 

4.23 

17.97 

81.51 

0.2500 

0.738 

HIRES 

1853.86184 

90.96 

3.67 

-6.05 

13.12 

0.1870 

0.639 

HIRES 

1853.88098 



67.06 

32.55 

0.1230 

0.646 

HIRES 

1879.90065 

-62.91 

4.21 

-11.44 

32.23 

0.1220 

0.095 

HIRES 

1880.99970 

1.46 

3.24 

10.18 

16.84 

0.1500 

0.494 

HIRES 

1882.08871 

76.62 

7.72 

-77.06 

46.76 

0.4450 

0.890 

HIRES 

1972.84557 

102.30 

4.25 




0.849 

HIRES 

2139.06062 

-134.45 

3.86 

-0.66 

15.61 

0.1360 

0.212 

HIRES 

HAT-P-53 

1877.91488 

217.41 

4.13 

3.37 

8.66 

0.1820 

0.708 

HIRES 

1879.81971 

220.06 

3.91 

-8.64 

5.38 

0.1710 

0.679 

HIRES 

1879.83898 



-6.30 

5.05 

0.1700 

0.689 

HIRES 

1880.95723 

-228.59 

3.82 

4.64 

4.22 

0.1740 

0.259 

HIRES 

1944.84748 

187.44 

5.04 

11.55 

4.67 

0.1010 

0.829 

HIRES 

1960.79667 

43.84 

4.68 

-10.07 

13.55 

0.1430 

0.959 

HIRES 

1972.75173 

-68.38 

4.23 

-14.42 

7.85 

0.1570 

0.054 

HIRES 


Note. — Note that for the iodine-free template exposures we do not measure the RV but do measure the BS and S index. Such template 
exposures can be distinguished by the missing RV value. 

^ The zero-point of these velocities is arbitrary. An overall offset 7rei fitted to these velocities in Section l3.3l has not been subtracted. 

^ Internal errors excluding the component of astrophysical jitter considered in Section |3.3l 
^ Chromospheric activity index calculated following llsaacson fc Fischeil 1 I 20 T 0 II . 












